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We study neutrino process nucleosynthesis in massive stars using newly calculated cross sections, 
an expanded reaction network, and complete and self-consistent models of the progenitor star. We 
reevaluate the production of light isotopes from abundant progenitors as well as that of rare, heavy, 
proton-rich isotopes. In particular, new results are given for "B, i^F, i3«La, and i«°Ta. The 
production of these isotopes places limits on neutrino spectrum and oscialltions. 
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Nuclei can be synthesized in the mantle of a core col- 
lapse supernova by the neutrino process - energetic 
supernova neutrinos excite nuclei above particle breakup 
through neutral- and charge-current reactions, creating 
new daughter nuclei. While typically only 1 % of man- 
tle nuclei experience inelastic neutrino reactions, certain 
rare nuclei, one mass unit below abundant parent nuclei 
such as C, O, and Ne, may nevertheless be produced 
dominantly by the neutrino process. The final abun- 
dance depends not only on the instantaneous yield of the 
daughter isotope, but also whether that isotope survives 
subsequent processing. In many cases the heating asso- 
ciated with the passage of the shock destroys the daugh- 
ter isotope: the surviving nuclei may be only those pro- 
duced post-shock, after the relevant shell has expanded 
and cooled. 

Two recent developments make a re-examination of the 
i/-process timely. First, we now know about neutrino 
oscillations, which could alter certain j/-process yields 
by enhancing charged-current production channels. Sec- 
ond, new data on the abundances of B and F - two key 
z/-process products - have been obtained. Prochaska. 
Howk. and Wolfe recently observed over 25 elements in 
a galaxy at redshift z — 2.626, whose young age and high 
metallicity implies a nucleosynthetic pattern dominated 
by short-lived, massive stars. Their observation of a so- 
lar B/0 ratio in an approximately 1/3-solar- metallicity 
gas argues for a primary (metal-independent) production 
mechanism, such as the i^-process (making ^^B), rather 
than a secondary process, such as cosmic ray proton spal- 
lation reactions on interstellar CNO seed nuclei (making 
^"B and ^^B). The new F abundance data of Cunha et al. 
showing a low F/0 ratio in two lo Centauri stars argue 
against AGB-star production of F. and are quite consis- 
tent with i^-process models (though also with production 
in cores of stars sufficiently massive to be Wolf-Rayet 



stars at the beginning of He burning) . 

Until the present effort the most complete neutrino 
process calculations were those of Q , who evaluated pro- 
ductions in a 20 Mq Pop I star evolved without mass loss, 
including semiconvection, and using the Caughlan et al. 
i2c(a, 7)1^0 reaction rate. The nuclear chemistry of 
the coproduccd protons and neutrons, which proved to 
reduce productions of important isotopes like ^-'^B, ^^N 
and ^^F, was followed in a nuclear reaction network. 
The effects of shock-wave heating and post-shock neu- 
trino process production, in shells expanding off the star 
and cooling, were evaluated. Later Timmes et al 
extended this calculation to a full galactic model, inte- 
grating the neutrino process over a range of progenitor 
stars with evolving metalicity. 

This letter extends this earlier work in several impor- 
tant ways. One is the incorporation of mass loss in the 
evolution of the progenitor star. Second, for the first 
time a reaction network is employed that includes all of 
the heavy elements through bismuth using updated reac- 
tion rates 0: the work of 0,0] ranged only up to zinc. 
This allows us to add selected neutrino reactions - the in- 
clusion of neutrino cross sections for the entire extended 
network has not yet been attempted - for products that 
may serve as electron- neutrino "thermometers." These 
cross sections are evaluated in a model that we believe is 
suitable for heavy nuclei. Third, the nuclear evaporation 
process - emission of a proton, neutron, or a ~ is treated 
in a more sophisticated statistical model that takes into 
account known nuclear levels and their spins and parities. 

The various partial neutrino cross sections are calcu- 
lated as two-step processes, as in The charged- and 
neutral-current cross sections are evaluated as function 
of excitation energy in the final nucleus, then a statis- 
tical model is used to evaluate the subsequent decay by 
particle or 7 emission. For the p- and sd-shell nuclei 
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the Gamow- Teller (GT) contributions to the {ve, e^) and 
{i^, v') responses were taken from Qhw shell model diag- 
onalizations, as appropriate. The Cohen-Kurath (^''N) 
and Brown- Wildcnthal interactions were used. For 
^^C we adopted energy and GT strength of the main GT 
transition, to the T = 1 state at 15.11 MeV in ^^C or 
its analogue in ^^N. from experiment. The double-magic 
nucleus ^^O has no GT response in the Qhw limit. All 
other contributions to the neutrino cross sections have 
been determined within the random phase approxima- 
tion (RPA) considering multipoles up to J = 4 and both 
parities. The RPA model, described in 0, treats pro- 
ton and neutron degrees of freedom separately and em- 
ploys a partial occupancy formalism for non-closed-shell 
nuclei. The residual interaction is a zero-range Migdal 
force. As realistic shell-model calculations for the heavy 
nuclei (^^^^Ba, ^^^La, ^^^Ta, i*°Hf) are yet not practical, 
the entire response was calculated using RPA. 

In the second step we use the statistical model 
SMOKER [nj to calculate, for each final state with well- 
defined energy, angular momentum, and parity, branch- 
ing ratios for p, n, a and 7 emission. SMOKER uses ex- 
perimentally determined levels in the daughter nucleus, 
supplemented at higher energies by an appropriate level 
density formula [ll|. If the decay leads to an excited 
level of the daughter nucleus above particle threshold, 
the subsequent decay of this level is treated similarly. 
The yield of a given nucleus is obtained by folding the 
various branching ratios, as a function of energy, with 
the neutrino response function. 

This is qualitatively the same procedure used in 0. 
There full multi-shell shell-model calculations were done 
through -'^^O; in the sd shell, however, the first-forbidden 
response was taken from the simpler Goldhaber- Teller 
model. The sd positive-parity shell model calculations 
for ^''Mg, ^®Si, and '^^S were also truncated. Perhaps 
more important, the branching ratios were evaluated 
in with a statistical model that lacked the capabili- 
ties of SMOKER (e.g., experimental level densities and 
spin/parity selection rules). In general, the present cross 
sections turn out to be slightly smaller than those of 
0. Detailed partial cross sections for the heavier nu- 
clei (i38Ba, i39Ta, i^^Hf, ^^^Ta) have not been previously 
calculated. 

For the neutrino spectra we took Fermi-Dirac distri- 
butions with zero chemical potential and temperatures 
T = 6 MeV for ^ and r neutrinos and their antiparti- 
cles and T = 4 MeV for and Vf.. Very recent super- 
nova simulations [l^ find somewhat harder v^^r spec- 
tra with T = 5.9 MeV and a degeneracy parameter 
a = 1.5, predicting slightly larger average neutrino ener- 
gies and increased cross sections (for example by 13 % for 
^°Ne(z^,i//)). The Pg energies are also found to be more 
energetic than what we assumed, but they have only little 
influence on the nuclei studied here. 

Potentially important i^-process candidates can be 



FIG. 1: Dominant neutrino process cross sections for produc- 
tion for ^^B, '^^F, '^^^La, and ^^'^Ta as a function of neutrino 
temperature (degeneracy parameter a — 0). 
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FIG. 2: Electron neutrino charged current cross sections on 
^^*Ba for cascades up to two particle emission. 



^ 2 

CM 

E 

O 

CM 
I 

O 

\ 
c 
_o 

u 

0) 

m 

in 
m 
o 
u 




-4 

4 6 8 

neutrino temperature / MeV 



identified by looking at the abundances in the star dur- 
ing the explosion, after the passage of the shock. This 
includes in particular radioactive parent nuclei, existing 
in nuclear statistical equilibrium or produced in the pas- 
sage of the shock wave. Table m shows some of the can- 
didate reactions that were identified from abundances 
found in our progenitor star 3.85 s after core bounce, 
the time when the shock reaches the base of the he- 
lium shell, located at radius 4 x 10^ cm and mass co- 
ordinate 6.3 Mq. The table includes an estimate of the 





TABLE I: Heavy i^-process candidate reactions as derived 
form a 25 stellar model 3.85 s after core bounce. 0-0,42 
is the cross section in lO"'''^ cm^ that would be required for 
solar production of the isotope. 



product 


CQ,42 


p circnt 


process 




51y 


3.0 


52Fe 


(i^e,e" n) 




^^^Mn 


2.0 


56Ni 


(!Ve,e" p) 




^»Kr 


11; 29 


»OKr 


iy^v' 2n); (;^e,e~ 


2n) 


«*Sr 


19; 50 


8«Sr 


(i/,//'2n); {uE,e~ 


2n) 




102 


85Rb 


{ue,e~ n) 




"8Ce 


33; 88 


"°Ce 




2n) 


i»OTa 


58 


182 jjf 


{ve,e~ 2n) 






118 


ISljjf 


(i/e,e~ n) 






22 




(z7g,e+ 2n) 




196 jjg 


74 


198 


(1^,1^' 2n) 






191 


198 


(i^e,e~ 2n), (!7e,e 


+ 2n), ... 



neutrino cross section aQ, required to produce the so- 
lar abundance of each isotope Yq, measured with re- 
spect to ^^O, when averaged over all ejecta: ctq := 
47ry0A^cjccta/?^ly /^j^^ja^ i^parontr"^dm. In this integration 
we use the progenitor radial coordinates r to locate the 
various shells, which would underestimate the needed 
cross sections for isotopes produced dominantly after 
shock wave passage (when the shell is expanding out- 
ward). We note the two most favorable candidates 
(smallest required cross sections), ^^V and ^^Mn, were 
produced at significant levels in the i^-process work of 
01 : those calculations used an interpolation method to 
generate neutrino cross section estimates for all nuclei in 
the network. (However, the ^^V and ^^Mn were made as 
^^Mn and ^^Co by the neutral current reactions {ly, ly'p) 
off of ^^Fe and ^^Ni, respectively, not by the charge- 
current reactions listed in the table.) The remaining re- 
actions in Tabic will be explored in future, once the 
needed neutrino cross sections are evaluated. 

Our full i^-process calculations included all dynamics, 
explosive and 7-process nucleosynthesis, and all net- 
work reactions destroying parent or daughter isotopes 
involved in i/-proccss synthesis. Supernova calculations 
were done with the KEPLER code [l^ , starting from the 
15 and 25 M© models S15 and S25 of 0. We assume a 
total energy of 5 x 10^^ erg per species with a luminos- 
ity exponentially decaying after onset of core collapse on 
a time-scale of 3 s and constant neutrino temperature of 
4 MeV for the i/^ and and of 6 McV for the z/^, Df^, I'r, 
and Dt, with zero degeneracy parameter (a = 0). 

We evaluated neutral current reactions and two par- 
ticle cascades for ^^c, ^^N, i^.isq^ 20,22]s^g^ ^^'^^^Mg, 
28,30gj^i39]^g^ and ^'^^Ta. For charged current reactions we 
include two particle cascades for ^^^Ba and ^^"Hf (Fig- 
ure Figure El shows that the 2n channel is important 
for these nuclei. 

It was recently argued by Goriely et al.^^, foUow- 
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FIG. 3: Production of "*La in a 25 M© star. We give the 
production of ^"^^La without neutrinos, with the charged cur- 
rent reaction ^■^*Ba(!-'c,e~ ), with the neutral current reaction 
^''^La(i/,i/'n) and with both reactions. Additionally, we show 
the SN shock temperature, the pre-SN density, and the neu- 
trino energy fluence Fe,v assuming a total neutrino energy of 
3x10^5^ erg. 
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TABLE II: Production factor relative to solar normalized to 
^^0 production as a function of T^. (for charged current only) 
and using 6 MeV for the /i and r neutrinos. 



star 


product 


(no i^) 


(no Uc) 


4 MeV 


6 MeV 


8 MeV 




"B 


0.011 


1.509 


1.899 


3.291 






15 N 


0.396 


0.480 


0.486 


0.530 




15 Mq 


19 p 


0.375 


0.577 


0.643 


0.914 






i^»La 


0.190 


0.279 


0.974 


1.734 


2.456 




i^opa 


0.599 


1.016 


2.751 


4.628 


6.026 




"B 


0.004 


0.828 


1.170 


2.384 






15 N 


0.039 


0.112 


0.118 


0.157 




25 M© 


19 p 


0.105 


0.300 


0.366 


0.643 






"«La 


0.106 


0.192 


0.901 


1.604 


2.244 




i«»Ta 


1.382 


2.360 


4.238 


6.238 


7.102 



ing the suggestion first made in that charged cur- 
rent reactions on ^'^^Ba would dominate the production 
of -'^^^La. We confirm this, finding the 7-process contri- 
bution is small and the neutral current reaction on ^'^^La 
is insignificant (Figure O. Sufficient ^'^^La is made for 
about solar co-production with ^^O in massive stars (Ta- 
blcini). The key is the enhancement of ^'^^Ba by the 
s-process prior to the SN (Figure^. Conversely, charge- 
current production of ^^^Pt is insignificant because of the 
destruction of the parent ^^^Ir in the s-process. 
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FIG. 4: Production of La in a 25 M© star and its neutrino 
process progenitor nuclei, ^^^La (neutral current) and ^^^Ba 
(charged current). The mass fraction of these isotopes as a 
function of the enclosed mass is shown before (gray) and after 
(black) the supernova explosion. 
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TABLE IIL Production factor relative to solar normalized 
to ^^0 production as a function of fj, and r neutrino tem- 
perature (neutral current) and using 4MeV for the electron 
(ani-)neutrinos (for charged current only). "Hax" are the re- 
sults from 3 using Haxton's u cross sections and "Kol" for 
the new rates of this paper by Kolbe. 

15 Mq 25 Mq 

product _ 6 MeV _ _ 8 MeV _ _ 6 MeV _ _ 8 MeV _ 





Hax 


Kol 


Hax 


Kol 


Hax 


Kol 


Hax 


Kol 


"B 


1.65 


1.90 


3.26 


3.99 


0.95 


1.17 


1.36 


1.85 


19 p 


0.83 


0.64 


1.28 


0.82 


0.56 


0.37 


1.03 


0.55 




0.46 


0.49 


0.54 


0.58 


0.09 


0.12 


0.15 


0.19 


i38La 




0.97 




1.10 




0.90 




1.03 


i«OTa 




2.75 




3.07 




4.24 




5.24 



apparent overproduction. 

The ^^B yield is somewhat higher than i (Table HH), 
but we note that this may also vary by up to a factor 2 
when modifying the ^^C(q;,7) rate within its current 30 % 
range of uncer tainty 0. ^^F is typically 50 % lower (Ta- 
blcnj due to the reduced cross section. The neutral cur- 
rent cross section is uncertain as it depends sensitively on 
the strong quenching of the GT strength as predicted by 
the shell model and what fraction of the GT strength re- 
sides above the particle threshold; experimental guidance 
is needed. The neutrino contribution to ^^N is increased 
by about 50% using our new rates, but the total yield 
remains low. 

The ^^^La and ^®°Ta yields are very sensitive to the 
z^e temperature, with ^'^^La the better thermometer as 
it is made almost exclusively in this channel (Tables ^ 
and IIII|) . Such a thermometer is potentially quite im- 
portant because of the role of the as yet unmeasured 
neutrino mixing angle 813. This parameter, crucial to 
proposed terrestrial long-baseline neutrino searches for 
CP violation, will generate matter-enhanced ^ os- 
cillations at a density corresponding to the atmospheric 
6m} . Naively this corresponds roughly to the base of the 
carbon zone. Such an oscillation would in turn "heat" the 
fgS (assuming a normal rather than inverted mass hier- 
archy), thus enhancing charged current rates. As most 
of surviving "'^^^La is produced at higher densities, in the 
outer half of the neon zone, one might conclude that its 
production does not probe this interesting physics. How- 
ever, it has been recently noted that v — v scattering may 
push the MSW resonance toward higher densities, thus 
into the neon zone jl7| . This is clearly an interesting 
question deserving further study. 
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The 7-process and neutral-current I'-process each ac- 
count for about 25 % of ^^''Ta, with the remainder coming 
from the charged current v process. It is somewhat over- 
produced. However, our calculations do not distinguish 
production in the 9~ excited state from production in 
the l"*" ground state, while only the former is long lived. 
Estimates indicated thermal freeze-out should leave 
30-50 % in the isomeric state. A significant part of what 
is made by neutrino process from low-spin parent nu- 
clei, however, might preferentially cascade down to the 
l'*" state as the v interaction is dominated by low multi- 
poles (A. Hayes, priv. com.). This could account for the 
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